Mon. Not. R. Astron. Soc. 0Q0.mfT8l(2009) Printed 23 July 2009 (MN BTgX style file v2.2) 



Galaxy Zoo Green Peas: Discovery of A Class of Compact 
Extremely Star-Forming Galaxies * 

Carolin Cardamone 1 ' 2 !, Kevin Schawinski 2 ' 3 , Marc Sarzi 4 , Steven P. Bamford 5 , 
Nicola Bennert 6 , C. M. Urry 2 ' 3 , Chris Lintott 7 , William C. Keel 8 , John Parejko 9 , 
^ ; Robert C. Nichol 10 , Daniel Thomas 10 , Dan Andreescu 11 , Phil Murray 12 , M. Jor- 

; dan Raddick 13 , Anze Slosar 14 , Alex Szalay 13 , Jan VandenBerg 13 

, 1 Astronomy Department, Yale University 208121, New Haven, CT 06520, U.S.A. "~ 

i 2 Yale Center for Astronomy and Astrophysics, Departments of Physics and Astronomy, Yale University, New Haven, CT 06520, USA 

^ ■ 3 Department of Physics, Yale University, P.O. Box 208121, New Haven, CT 06520, USA. 

5 ' 4 Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts ALIO 9AB, UK. 

■ 5 'Centre for Astronomy and Particle Tlieory, University of Nottingham, University Park, Nottingham, NG7 2RD, UK. 

■ ^Department of Physics, University of California, Santa Barbara, CA 93106, USA. 
7 Department of Physics, University of Oxford, Oxford 0X1 3RH, UK. 

1 | 8 Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL, 35487, USA. 
, 9 Department of Physics, Drexel University, Philadelphia, PA 19104, USA. 

U 1( 'institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth, POl 2EG, UK. 
"LinkLab, 4506 Graystone Ave., Bronx, NY 10471, USA. 

(-H ' 12 Fingerprint Digital Media, 9 Victoria Close, Newtownards, Co. Down, Northern Ireland, BT23 7GY, UK. 

f*) 13 Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA. 

I i 14 Berkelev Center for Cosmological Physics, Lawrence Berkeley National Lab, 1 Cyclotron Road, MS 50-5005, Berkeley, CA 94720, USA 

o 

CO 

^ 1 23 July 2009 

ABSTRACT 

I/"*) , We investigate a class of rapidly growing emission line galaxies, known as "Green Peas," 

■ first noted by volunteers in the Galaxy Zoo project because of their peculiar bright green 

colour and small size, unresolved in SDSS imaging. Their appearance is due to very strong 
optical emission lines, namely [O III] A5007 A, with an unusually large equivalent width 

r***> ' of up to ^1000 A. We discuss a well-defined sample of 251 colour-selected objects, most 

of which are strongly star forming, although there are some AGN interlopers including 8 
newly discovered Narrow Line Seyfert 1 galaxies. The star-forming Peas are low mass galax- 
ies (M~ lO 8 ' 5 - 10 10 M Q ) with high star formation rates (~ 10 M Q yr _1 ), low metallicities 
(log[0/H] + 12 ~8.7) and low reddening (E(B - V) ^ 0.25) and they reside in low density 
environments. They have some of the highest specific star formation rates (up to ~ 10~ 8 
r>J ' yr _1 ) seen in the local Universe, yielding doubling times for their stellar mass of hundreds of 

Myrs. The few star-forming Peas with HST imaging appear to have several clumps of bright 
star-forming regions and low surface density features that may indicate recent or ongoing 
mergers. The Peas are similar in size, mass, luminosity and metallicity to Luminous Blue 
Compact Galaxies. They are also similar to high redshift UV-luminous galaxies, e.g., Lyman- 
break galaxies and Lyman-a emitters, and therefore provide a local laboratory with which to 
study the extreme star formation processes that occur in high-redshift galaxies. Studying star- 
bursting galaxies as a function of redshift is essential to understanding the build up of stellar 
mass in the Universe. 

Key words: galaxies: evolution, galaxies: formation, galaxies: starburst, galaxies: dwarf, 
galaxies: high-redshift, galaxies: Seyfert 
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1 INTRODUCTION 



2 DATA 



In this paper, we report on the discovery of an intriguing class of ob- 
jects discovered by the Galaxy Zoo project. The Galaxy Zoo project 
dLintott et aT1l2008h has enlisted the help of over 200,000 mem- 
bers of the public to morphologically classify almost 10 galax - 
ies from the Sloan Digital Sky Survey (SPSS; lYork et al .1 2000). 
The Galaxy Zoo websitj]] provides a randomly selected gri com- 
posite colour image from the SDSS main galaxy sample and asks 
the volunteers to classify the morphology of the displayed ob- 
ject. One advantage of this distributed approach to classification 
is the fact that each object will receive multiple, independent clas- 
sifications, and so one can treat the distribution of classifications 
for each object in a statistical sense. These classifications have 
led to numerous results in galaxy formation and cosmology (e.g 



'logy (. . 

Skibbaetal. 2008; Land et al. 2008; Bamford et al. 2009; Lintott 



20091 ; ISchawinski et al.ll2009l : ISlosar et al.ll2009l ; lDarg et alj 2009). 

In addition to the website used for classification, Galaxy Zoo 
also provides an online discussion forun{3 where volunteers may 
ask questions about unusual or challenging objects. This allows us 
to tap into another advantage of human classifiers: they can easily 
identify and then investigate odd objects. One such class of highly 
unusual objects was named 'Green Peas' as they appeared to be 
unresolved round point sources that looked green in the gri com- 
posite. FigureQjshows three example Pea images found by users, as 
well as a more typical galaxy at the same redshift (z~0.2). The vol- 
unteers rapidly assembled over a hundred of these objects in a ded- 
icated discussion threacjfl Most of these were class ified as stars in 
the SDSS photometric pipeline dLupton et alj|200lh . It quickly be- 
came apparent that these objects represent a distinct group: they all 
had galaxy-type spectral features (as opposed to broad-line quasar 
spectra or stellar spectra), and their green colour was driven by a 
very powerful [OIII] A5007 Aemission line that substantially in- 
creased the r-band luminosity relative to the adjacent g and i band 
(green being the colour represented by the r-band in the SDSS 
colour composites). As a result of this selection, the Peas are found 
at redshifts 0.112 < z < 0.360, mostly beyond the main galaxy 
sample but much nearer than luminous quasars. This discovery 
prompted our investigation into the nature of these small [O III]- 
emitters. 

In Section [2] we present the sample selection, and analyse 
their photometric properties (§ 12. U . space density (§ 12. 2t . mor- 
phologies in SDSS (§ |2,3I >, environments (§ |2.4| l, spectral properties 
( i]2. 51 - ^2. 6b and HST imaging morphologies (§ 12.71 1, The Peas turn 
out to be largely star-forming objects with some AGN interlopers. 
We look into the nature of the few AGN Peas in Section[3] We study 
the nature of the star-forming galaxies in Section [4] and compare 
them to other known samples of galaxies in Sectionf5] Throughout 
this paper we assume Ho = 71 kms -1 Mpc -1 , fi m = 0.3 and 
JIa = 0.7, consistent with the Wilkinson Microwave Anisotropy 
Probe 3-year results in combination with other cosmology probes 
dSpergel et alj2007h . 



http://www.galaxyzoo.org 

2 http://www.galaxyzooforum.org 

3 We wish to thank the "Peas Corps" for "giving Peas a chance:" includ- 
ing, Elisabeth Baeten, Gemma Coughlin, Dan Goldstein, Brian Legg, Mark 
McCallum, Christian Manteuffel, Richard Nowell, Richard Proctor, Alice 
Sheppard, Hanny van Arkel. 



2.1 Sample Selection 

Our sa mple of Peas is taken from the SDSS DR7 spectroscopic 
sample dAbazaiian et~afll2009r) . A survey of a quarter of the sky, 
the SDSS provides phot ometry of 357 million unique objects in five 
filters, it, g, r, i and z jFukugita et al. 1996) and spectroscopy of 
many objects. Using the CasJob application provided by SDSS, 
we uniformly searched the DR7 spectroscopic sample for Peas 
(originally noticed by eye in Galaxy Zoo) in the redshift range 
0.112 < z < 0.360 where the [O III] A5007 Aline is in the r- 
band filter. 

To define colour selection criteria, we compared the sample 
of ~100 Peas identified by the Galaxy Zoo volunteers to a com- 
parison sample of 10,000 galaxies and 9,500 QSOs at the same 
redshifts over the colour space defined by the 5 SDSS bands. The 
10,000 galaxies were selected to match the redshift and g-band- 
magnitude distributions of the Peas. The QSO sample contains all 
spectroscopically confirmed QSOs in the Peas' redshift range, be- 
cause the QSOs are overall too luminous to match the Peas mag- 
nitude distribution. Figure [2] displays two colour-colour plots with 
Peas (green crosses), comparison galaxies (red points) and compar- 
ison QSOs (purple stars). The green Pea colour selection is shown 
by the darkened black lines. The precise colour cuts were selected 
to avoid both the QSO and overall galaxy sequences and to high- 
light the unusual objects selected by eye in the Galaxy Zoo forum. 
The colour limits are: 



r - i < -0.2 
r - z sC 0.5 
g — r ^ r — 0.5 
u — r ^ 2.5(r — z) 



(1) 

(2) 
(3) 
(4) 
(5) 



We illustrate the effectiveness of this colour selection in Figure 
[2] It divides the Peas from the loci of both the galaxy and the 
quasar populations. This colour selection technique effectively uses 
the narrow-band survey technique common in high redshift galaxy 
searches, only here applying it to the broader SDSS filters. Because 
we are using broad filters, we are only sensitive to galaxies with ex- 
treme [O III] equivalent widths. Further, because we are using the 
r-band filter for colour selection, we are not sensitive to Pea-like 
objects at lower and higher redshifts. As seen in Figuref2] the Peas 
do indeed have distinct colours from the SDSS main galaxy sam- 
ple. This is especially noticeable in the r — i and g — r colours used 
to create the 3-colour images. Using the colour selection, we find 
a sample of 251 Peas taken from all SDSS spectroscopic galaxies 
whose spectra we further analyse to understand their properties. 

Although most of the Peas were identified by the SDSS spec- 
troscopic pipeline as having galaxy spectra (4 were classified as un- 
known), only 7 were targeted by the SDSS spectral fibre allocation 
as galaxies. Most were targeted as serendipitous objects, with the 
majority flagged as SERENDIP_BLUE, SERENDIPJ3ISTANT 
and TARGE T_QS0_FAI NT. These target flags are for objects lying 
outside the stellar locus in colour space and DISTANT here refers 
to distance from the stellar locus, QSCLFAINT is also used for ob- 
jects flagged as stellar that are both f ainter than i=20 mag a nd out- 
side the stellar locus in colour space tStoughton et al. 2002). These 



http : / /cas jobs . sdss . org/ Cas Jobs/ 
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Figure 1. Example g,r,i composite colour 50"x 50" SDSS images classified by Pea hunters, with the r-band representing green light. The distinctly green 
colour and compact morphology makes the Peas (left 3 images) easily distinguishable from the classical elliptical (right image). The elliptical galaxy is clearly 
red and has a smooth profile, while the Peas are r-band dominated and unresolved in these images, appearing like stellar point sources. All objects shown here 
are at z~0.2. 




Figure 2. Left: r — ivs g — r colour-colour diagram for 251 Peas (green crosses) and a sample of normal galaxies (red points) matched in z and g-band 
magnitude, and all QSOs (purple points) which lie in the same redshift range as the Peas, 0.112 < z < 0.36. Right: r — z vs u — r colour for the same 
classes. The Peas are most cleanly separated in the r — i, g — r colour diagram, where they stand out as particularly bright in the r-band. The colour cuts were 
selected to include the Peas identified by the Galaxy Zoo volunteers and to exclude both galaxies and QSOs. 



objects were targeted by fibers as they became available in a given 
field, so their selection function is not well determined. Without a 
uniform selection for the Peas across the sky, their absolute space 
density cannot be accurately assessed. 



2.2 How common are the Peas? 

Because the spectroscopic selection is biased in unquantifiable 
ways, i.e., based on the availability of a spectral fibre in a given 
pointing, the space density of the Peas is difficult to assess. In or- 
der to estimate the space density of the Peas in SDSS, we need 
to search the entire SDSS photometric catalogue for objects with 
our colour selection criteria. We note that the Pea colour selection 
can also return much higher redshift objects by finding other emis- 
sion lines in the r-band filter. Therefore we first search the spectro- 



scopic sample to help understand these contaminants. Dropping the 
Pea redshift selection limits, we find 198 objects at higher redshift 
in the spectoscopic sample which fall into our colour-selection re- 
gion. These are mostly QSOs (only 4 have a spectral type labelled 
as galaxy by the SDSS pipeline), which cluster at redshifts z=1.2, 
where the 2800A Mg II line falls into the r-band, and z=3.0, where 
the 1546 A CIV line falls into the r-band. Very roughly, in the en- 
tire spectroscopic database, there are comparable numbers of Peas 
and higher-z QSOs in the colour-selection region. Searching the en- 
tire SDSS photometric catalogue (PhotoOb j) regardless of spec- 
troscopic information, we limit ourselves to objects with the same 
r-band magnitude range as the Peas (18 ^ r ^20.5) and to ob- 
jects with similar compactness (petroracLr ^ 2.0 "), in order to 
limit the contaminants from other redshifts. We further added the 
requirement of small g, r and i band photometric errors to avoid the 
scattering of contaminants with poor photometry into the colour se- 
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lection region. This search returns 40,222 objects. The unique area 
of the SDSS DR7 footprint covers 8423 sq. deg, so this implies 
a rough spatial density estimate of 5 per square degree. Strictly 
speaking this is an upper limit because our selection likely still con- 
tains QSOs from higher redshifts. Looking at QSO number counts, 
we would expect to s ee ~3 per sq. degree in our magnitude range 
dRichards et akOoOsh . leaving two Peas per square degree brighter 
than 20.5 magnitudes. Therefore, we conclude the Peas are rare ob- 
jects. 



2.3 Morphology from SDSS imaging 

Of the 251 Peas in our sample, 215 were classified as morpholog- 
ical type STAR and not as extended objects by the SDSS pipeline. 
Compared to the size of the galaxies as measured by the SDS S 
Pipeline (Petrosian radius: iBlanton et al.l200ll ; lYasuda et alj20oih . 
they are significantly smaller. This is expected because the spectro- 
scopic galaxy sample was selected to be both resolved and brighter 
than the faintest Pea dStrauss et ai]|2002l) . The typical resolution of 
the SDSS images is large (J> 1"), just below the peak of the sizes 
of the Peas as measured by the SDSS Petrosian radius. At the typi- 
cal redshifts of the Peas, this angular scale corresponds to an upper 
limit on the physical half-light radius of approximately 5 kpc. 



2.4 Environment 

Because the Peas are at redshifts higher than the main spectroscopic 
galaxy sample, we measure the projected densities around the Peas 
and a sample of random galaxies matched in both luminosity and 
redshift, counting the number of projected neighbours within a ra- 
dius of 1 Mpc and brighter than an absolute i-band magnitude of 
-20.5 at the redshift of the target are counted. This magnitude limit 
was selected to be well above the detection limit of the highest red- 
shift Peas. Figure[3]shows the cumulative distribution of neighbour 
counts for the Peas (solid line) and the control sample of galaxies 
(dotted line). Although this is a simple test, neglecting foreground 
and background contamination in the neighbour counts, and hence 
underestimating the difference between high and low density envi- 
ronments, it is clear that the Peas inhabit significantly lower den- 
sity regions than typical galaxies of the same i-band luminosity. 
A Kolmogorov-Smirnov test indicates that the difference is signifi- 
cant at a greater than 8-cr level. We find the median environmental 
density around the Peas is less than two-thirds of that around nor- 
mal galaxies. 



2.5 Spectral Analysis 

We downloaded all of the Peas' spectra from the DR7 archive 
and carefully re-fit them, paying close attention to both the con- 
tinuum and emission lines. We fit the stell ar continuum an d 
ionised-gas emission, following the technique o flSarzi et al. I fcood) . 
using the corresponding PP XF (PENALISED PIXEL FITTING; 
ICappellari & Emselleml [20041) and GANDALF (GAS AND AB- 
SORPTION LINE FITTING) IDL (Interactive Data Language) codes 
adapted for dealing w ith SDSS datfl Stellar population templates 
jTremonti et"afl f2004) and Gaussian emission-line templates were 
simultaneously fitted to the data. When it improved the overall fit (9 




150 



N(r < IMpc) 



Figure 3. The cumulative distribution of neighbour counts for the Peas 
(solid line) is significantly different than a comparison galaxy sample (dot- 
ted line). No Peas live in the highest density environments. 



cases), the Gaussian emission line templates included both broad- 
line and narrow-line components for the Balmer Series. 

To ensure acceptable fits, we limit our spectroscopic sample to 
those Peas with a S/N greater than 3 in the continuum near the Hf3 
and Ha regions (specifically, measuring the S/N in the bands 6350- 
6500A and 5100-5250A). Additionally, we limited our emission- 
line classification sample to those Peas with a S/N ^ 3 detection 
in each of the emission lines; Ha, H/3 [O III ] A5007 Aand [N 

II] A6583 A, following iKauffmann et all d2003l) . We note that for 
many of our objects near z~0.3, sky lines fall on top of the [O 

III] line and the Hf3 line, and we removed all of these objects from 
our sample. One of these objects has broad Balmer lines, with both 
sky lines falling inside the H/3 profile. This objec t is identified as a 
Narrow Line Seyfert 1 (NLS1) in the literature dZhou et alJboOfj) . 
but we do not consider it further in this paper. These cuts result 
in a sample of 103 narrow-line objects to be further analysed (see 
§ 12, 6t plus nine objects best fit by a two-component Gaussian in the 
Balmer lines. The width of the broad Gaussian components ranged 



5 Both codes can be downloaded 

http : / /www . strw . leidenuniv .nl/sauron/ 



from 



from just over 600 km s~ to ~5000 km s~ . Eight of these ob- 
jects have FWHM narrower than 2000 km s _1 and are thus classi- 
fied as Narrow Line Seyfert Is (see Sj3). These eight SDSS spectro- 
scopic objects have not previously bee n identified as NLSls, (e.g., 
Willi ams et al.l2002l ; lzhou et al.l2 006). Three examples of spectral 
fits are shown in Figure [4] The top 2 spectra are of narrow line ob- 
jects and the bottom spectra shows a NLS 1 object. In all three cases 
one can easily see a prominent [O III] doublet near 5007 A. 

In summary, we have 9 Peas with two component fits, (1 broad 
Line Seyfert 1 and 8 Narrow Line Seyfert Is) and 103 narrow line 
Peas. 



2.6 Spectral Classification 

The SDSS spectra cover the observed range 3800-9200 A at a re- 
solving power of R~1800. At the Peas' redshifts, this range in- 
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3500 4000 4500 5000 5500 6000 6500 

Rest Wavelength [A] 

Figure 4. Example spectral fits from GANDALF. The top two plots show 
typical star-forming Peas and the bottom plot shows a typical Narrow Line 
Seyfert 1. In black is the rest-frame observed spectrum and in red the fit 
from GANDALF. The SDSS filter band passes are included as blue dotted 
line, shifted into the rest-frame of the Pea. Notice in all examples, the [O 
III] A5007 A line is redshifted inside the r bandpass. The Peas show very 
strong emission lines, including clear lines detected throughout the Balmer 
series. 




4800 4900 5000 5100 



6400 6500 6600 6700 




Ha + [NII] region 




4800 4900 5000 5100 



6400 6500 6600 6700 
K, 



Figure 5. Example spectral fits from GANDALF in the regions around Ha 
(right) and H/3 (left). The top plots show a typical star-forming Pea and the 
bottom plots show a typical Narrow Line Seyfert 1 . 



Table 1. Spectroscopic Classification 



Type 



Number 



Broad Line AGN (Seyfert Type 1) 1 

Narrow Line Seyfert Is 8 

Narrow Line AGN (Seyfert Type 2) 10 

Transition Objects 1 3 

Star Forming 80 



Total 



112 



eludes the regions around both the H/3 and Ha spectral lines. Ex- 
ample fits to these lines are shown in Figure|5] 

We use a classical emission line diagnostic originally devised 
by Baldwin, Phi llips & Terlevich 1981 (hereafter BPT) and mod- 
ified by others {O sterbrock & Pogge 1985; Veilleux & Osterbrockl 
ll987l : lKewlev et alj|200ll, [2006; Kauff mann etai]|2003h to classify 
the narrow line objects. Emission line diagnostics probe the nature 
of the dominant ionising source and separate the galaxies dom- 
inated by ongoing star formation from those dominated by non- 
stellar processes (namely Seyfert and LINER galaxies). The classi- 
cal BPT diagram, which compares the ratio of [O III] A5007 A/H/9 
to [N II] A 6583 A/Ha, has been s hown to be an efficie nt measure 
of t he ionising source in a galaxy iKewlev et alj 2006). Addition- 
ally, iKewlev et all J200ll) calculated the maximum starburst con- 
tribution from theoretical models, including modern stellar popu- 
lation synthe sis, photoionisation and shock models (labelled KeOl 
in Figure [6]l. Kauffm ann et al.1 d2003h empirically shifted this line 
to separate purely star-forming objects from the rest using a set 



of ~ 23,000 SDSS spectra (labelled Ka03 in Figure |6j. Objects 
that lie in between the KeOl and Ka03 lines are transition galaxies 
containing a mix of star formation and a central AGN component 
(Kewley 2006). We note that the SDSS spectral fiber includes only 
the light of the central 3" of the galaxy, so emission originating 
from the central AGN-ionised regions can be mixed with emission 
originating from extended star formation. 

Analysis of the narrow-line Peas' spectra (Figure[6]l revealed 
that the majority of the objects are star forming (80, red stars), but 
there are also 10 Seyfert 2 (blue diamonds) and 13 transition objects 
(green crosses). The location of the star-forming galaxies in the top 
left of the BPT diagram indicates they likely have lower metallicity. 
This sample of galaxies is discussed in greater detail in §|4] In Table 
[T] we summarise the spectroscopic identifications. 
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log([NII] A 6583/Ha) 

Figure 6. The BPT diagram classifies 103 narrow-line Peas (all with S/N 
^3 in the emission lines) as 10 AGN (blue diamonds), 13 transition ob- 
jects ( green crosses) and 80 starbursts (red stars). Solid line: Kew lev et alj 
1200 ll) maximal starbu rst contribution (labelled KeOl); Dashed line: 
Kauffmann et al. 1 2003) line separating purely star-forming objects from 
AGN (labelled Ka03). The clustering of the starbursts in the top left corner 
of the plot indicates that they likely have low metallicity. 



2.7 Morphology from HST imaging 

Because the SDSS resolution is too low to measure the actual sizes 
of the Peas, we searched the HST archive, MASlffl finding 5 public 
images. The observations are summarized in Table|2] 

We reduced the ACS data, starting with the pipeline-reduced 
*flt* images and using MULTIDRIZZLE to remove cosmic rays and 
defects, correct for distortion, and improve sampling of the PSF 
with a final scale of 0.04"/pixel. In the case of the WFPC2 data, 
we started with the pipeline-reduced *c0f* images and used DRIZ- 
ZLE with parameters pix.scale=0.5 and pix.frac=0.8, leading to a 
final scale of 0.023 "/pixel for the PC chip and 0.05 "/pixel for the 
WF4 chip, respectively. Figure|7]shows all 5 public HST images of 
the Peas compared to a single typical ground-based SDSS image 
(bottom right). 

For the ground-based image (Fig|7J bottom right), the SDSS 
archive lists a Petrosian radius of l"8, roughly 8 kpc at the Pea's 
redshift (z=0.2832). This object is clearly unresolved when com- 
pared to the HST image of the same galaxy (top right) and the size 
as measured by the SDSS pipeline is therefore an upper limit. Of 
the 5 Peas with HST data, one is classified as a NLS1 (see ij3j bot- 
tom centre image), three are identified as star forming (see Sj4j top 
images) and the last one has sky lines over the [O III] region and 
is thus unclassified (see i]2.5l bottom left). The three star-forming 
galaxies were imaged as par t of a study of loc al UV-luminous 
galaxies (Heck man et alj200 5: Overzi er" et al.l2008h . All three Peas 
classified as starburst galaxies reveal complex structures too small 
to be resolved in ground-based imaging. The morphology of the top 
right object for example, shows several "knots" instead of one cen- 
tral component. These knots may be different star-forming regions, 
suggesting a morphology typical of merger events. Although the 

e Multimission Multimission Archive at STScI, 

http://archive.stsci.edu /index. html 



Peas live in low density environments, their star formation could 
still be driven by merging activity. The top centre image shows a 
central component with extended structures that look like two tidal 
tails reaching out to the east and south, possibly connecting to an- 
other galaxy in the south. The centre of the top left image seems 
also to consist of at least two knots rather than a smooth overall 
light distribution. The NLS1 (bo ttom centre) is from a study of 
AGN host galaxies dSchmitt 2006) and the last object (botto m left) 
was i maged serendipitously in a study of Kuiper belt objects jNollI 
l2007t) . In contrast to the star-forming galaxies, the NLS1, bottom 
centre, looks like a spiral seen edge on. This AGN aside, all of 
the HST images reveal complex structures much too small to be 
resolved in ground based imaging. Although our statistics are too 
low to make any conclusions on the general nature of the Peas, this 
is an interesting trend. 



3 THE PROPERTIES OF NARROW LINE SEYFERT 1 
PEAS 

NLSls constitute ~ 15% of low redshift (z<0.5) Seyferts 
dWilliamset ail [2002) and are characterised by H/3 line widths 
broader than classical narrow-line AGN but narrower than classi- 
cal Seyfert 1 galaxies. They appear to have Eddington ratios near 1 
and black hole masses below the typical Mbh-o" relation, suggest- 
ing a time delay between the growth of the gal axy and the growth of 
the central black hole mass (Ryan et al. 2007). Several authors use 
the [OIII] line as a surrogate for the bulge stellar velocity disper- 
sion a, assuming that the velocity field of the narro w-line region is 
domi nated by the stellar gravitational potential (e.g jBonning et all 
2005). However, the [OIII] line is known to often exhibit a blue 
wing that is ass ociated with gas out flow a nd this can affect mass 
determinations dMarconi et~al]|2008h . Thus. lKomossa et alj (2008) 
fit the blue wing separately and use only the main [OIII] component 
to estimate a. This careful fitting can decrease the galaxy host mass 
measu rement, placing the N arrow Line Seyfert 1 on the Mbh-c re- 
lation iKomossa et alj|2 008). 

Eight of the Peas are classified as NLSls, which we define 
as galaxies with 500 < FWHM H/S < 2000 km s" 1 . We de- 
termined the black hole and galaxy masses for our sample of 
N LSls. Black hole masses are estimated using the relation given 
bv lMcLure & Jarvisl J2002h : 

M BH = 10 7 ' 63 vl 000 LS 4 61 M @ , (6) 

where v 30 oo is the FWHM of the H/3 line divided by 3000 km/s 
and L44 is the luminosity at 5100A divided by 10 44 erg s _1 . To 
determine the galaxy mass, we fit the [O III] emission line with 
both a central narrow component and an additional blue wing. We 
then take the central narrow component as a proxy for the stellar 
velocity dispersion a. 

Figure [8] shows the stellar velocity dispersion of the bulge, 
estimated from uroni], vs black hole mass for the NLSls, nearly 
all of which lie well b elow the classical M-a relation (solid line, 
iTremaine et alj d2002l) : dashed lines are 1-a error contours) even 
though we fit the [O III] line with an additional blue wing com- 
ponent. This contrasts with a recent study that fo und this fitting 
metho d puts the NLS Is nearer to the M-<r relation dKomossa et alj 
2008). Instead, our results are consistent with studies m easuring the 
galax y mass using near-infrared bulge measurements dRvan et al.1 
2007). We also note that single broad-line AGN found in our sam- 
ple is consistent with the M-ct relation. Therefore the properties of 
our sample of NLSls are consistent with some of those found in 
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Table 2. HST Images 



SDSS Obj ID 


RA (J2000) 


Dec(J2000) 


z 


Instrument 


Chip 


Filter 


Exposure Time (s) 


Figure location 


587731187273892048 


351.41345 


0.75201 


0.2770 


WFPC2 


PC 


F606W 


3600 


top left 


588013384341913605 


141.50168 


44.46004 


0.1807 


ACS 


WF 


F850LP 


2274 


top centre 


587724199349387411 


10.22636 


15.56935 


0.2832 


WFPC2 


PC 


F606W 


3600 


top right 


587726879424118904 


344.49532 


-8.62438 


0.3081 


ACS 


WF 


clear 


1071 


bottom left 


587726032799400204 


211.91701 


2.29671 


0.3092 


WFPC2 


WF4 


F814W 


1200 


bottom centre 



SB 

587731187273892048 


SB 

588013384341913605 


* 


587726879424118904 


NLS1 — 



m -* ■ 

587726032799400204 


SB 

587724199349387411 

■ 



Figure 7. Hubble Space Telescope images of all Peas available in the archive. In each panel, we indicate a physical scale of 3 kpc at the redshift of the object 
(bar in top right of image). In the top row, we show three Peas classified as actively star-forming by their emission line ratios. In the bottom row, from left to 
right, an unclassified Pea, a Pea classified as a NLS1 and an SDSS image of a star-forming Pea at the same scale to highlight the degree to which the Peas are 
unresolved in typical SDSS imaging. These HST images illustrate that the upper limit on the Peas' physical scale based on ground-based SDSS imaging is a 
significant overestimate of the true physical size. Furthermore, these HST images reveal that the star-forming Peas exhibit a complex morphology indicative 
of significant disturbances that may be due to mergers and/or clumpy regions of star formation or extinction. 



the literature, and can help in the study of the location of NLS Is on 
the M-cr relation. 

To better characterise the nuclear emission using high energy 
data, we searched the database at HEASARC0 for additional data 
on our sample. Unfortunately, the Peas are well distributed through- 
out the 8,400 square degrees covered by the SDSS and not concen- 
trated in any of the areas covered by deep multi-wavelength sur- 
veys. None of the NLS Is are bright enough at soft X-ray wave- 
length s to be detected in the ROSAT All-Sky Survey dVoges et al.l 
1999). We note that the broad-line AGN is detected with a lumi- 
nosity of nearly 10 44 ergs s _1 , as is one of the type 2 Seyferts with 
an X-ray luminosity of nearly 3 x 10 44 ergs s _1 (0.2-2 keV). Given 
the redshift range of the Peas and the detection limits of the Rosat 
all sky survey, this limits the NLS Is' X-ray luminosity to below a 
few xl0 44 ergs s" 1 ; however, this upper limit is still well within 
the typical range of Seyfert luminosities. 

One of the HST images is of a NLS1 Pea (Figure [7] bottom 
centre). This galaxy looks distinctly different from the patchy HST 
images of the Peas powered by star formation. It appears to be an 
edge on disk, with no sign of morphological disturbance. This is 

7 NASA's High Energy Astrophysics Science Archive Research Center; 

http://heasarc.gsfc.nasa. gov/ 



consistent with what is seen in the morphologies of other samples 
of NLS Is tevan et alj20Q7l) . 



4 THE PROPERTIES OF STAR-FORMING PEAS 

From our spectral diagnostics we have 80 star-forming objects, 
with high signal-to-noise, which we look at in greater detail here. 
They are listed in Table [4] where we include information from the 
SDSS DR7 archive in Columns 1-5: SDSS Objld, RA, Dec, z and 
the equivalent width of [O III]. 



4.1 Star Formation Rates 

To determine accurate star-formation rates using the Ha fluxes, 
we corrected these recombination fluxes for both interestellar ex- 
tinction and fo r the underlying stellar absorption lines in the stel- 
lar continuum JKennicuttll 19981) . We measured the Balmer decre- 
ment, assuming an Ry = Ay/E(B — V) = 3.1, and using the 
ICardelli et all fl989) reddening curve and an intrinsic Ha / H/3 
value of 2.85 (the Balmer dec rement for ca se B recombination at 
T=10 4 K and n e = 10 4 cm~ 3 ; lLequeuxl2005h . There are a handful 
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Table 3. Narrow Line Seyfert 1 s 



SDSS Objld 


RA [degs] 


Dec [degs] 


z 


<7[OHI] [knis 1 ] a 


M BH [10 7 M Q ] b 


587726032799400204 


211.917006 


2.296711 


0.30920 


240.45 (4.47) 


0.100 (0.024) 


587733410983182549 


214.828750 


51.044473 


0.32363 


297.90 (8.69) 


0.429 (0.239) 


587733398646620415 


245.539241 


35.352080 


0.26601 


129.72 (2.47) 


1.270(0.323) 


587731521734640128 


117.387347 


28.568545 


0.33697 


226.91 (4.96) 


1.883 (0.599) 


587731892187037787 


172.279625 


57.934812 


0.31238 


305.33 (3.50) 


0.287(0.103) 


588017978351616137 


171.563561 


38.971510 


0.33651 


246.36 (4.38) 


0.102 (0.026) 


587739377230610665 


124.500818 


19.302802 


0.32452 


320.52 (3.83) 


1.081 (0.329) 


587739406266728813 


239.238254 


21.520959 


0.23314 


356.63 (3.31) 


0.720(0.139) 



a We measured the [OIII] line width fitting simultaneously for a blue wing and narrow Gaussian components. The measurements reported here are for the 

narrow Gaussian component with 3 sigma measurement error in parentheses. 
b The black hole masses are determined as in Eqn 1 , 1 sigma errors are in parentheses. 
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Figure 8. Using the [OIII] width as a measure of the host galaxy mass, we 
plot one broad-line AGN (star) and eight NLSls (diamonds) on the Mbh - 
a relation. The NLSls mostly lie far below the standard Mbh-o" relation, 
even though we determined galaxy mass from the narrow [OIII] component 
only. 



of star-forming Peas with Ha / H/3 less than 2.85, but these mea- 
surements are due to a combination of intrinsically low reddening 
and uncertainty in line flux determinations (Figure[9] top panel); we 
set the extinction equal to zero in these cases. Overall, we find that 
the reddening values for the Peas are low (Figure[9] bottom panel), 
with nearly all Peas having E(B — V) ^ 0.25. Finally, using the 
corrected Ha flu x measurements, we measured star-formation rates 
fcennicutd[l998h up to - 30 Meyr -1 . The measured star forma- 
tion rates are added to table|4]in Column 7. 
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Figure 9. The histograms of Ha / H/3 (top) and the colour excess, E(B — 
V), (bottom) as determined from the Balmer decrement. This distribution 
indicates that the line-emitting regions of star-forming Peas are not highly 
reddened, particularly compared to more typical star-forming or starburst 
galaxies. 



4.2 Stellar Mass 

The optical spectral energy distribution (SED) of the Peas is domi- 
nated by their strong nebular emission lines. Thus we cannot apply 
standard SED fitting techniques directly to their photometric data. 
Additionally, virial masses are inaccessible due to the low S/N of 
the spectral continuum and our inability to fit absorption lines. We 
therefore turn to the SDSS spectra, where the emission lines can 
be subtracted or blocked out where necessary. The continuum of 
the resulting emission-line-free spectrum has very poor signal-to- 
noise, so rather than fit this spectrum directly, we convolve it with a 



set of 19 medium-band filters jTaniguchill200i) and treat the result 
as medium-band photometric data. In other words, we construct an 
SED from the spectral continuum. To this, we add the GALEX near- 
and far-UV photometric data points where available (§|43) and fit 

the SED with stellar population models. 

W e employ a method similar to that of ISchawinski et al.l 
j2007l) and model the star formation history with two bursts, to ac- 
count for the possible presence of an underlying old stel lar popula- 
tion. We use the stellar models o f lMarastorj Jl 9981, 120051) with the 
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Salpeter IMF and a ran ge of metallicit i es. W e also account for dust 
extinction following the lCalzetti e t al. (2000) law and fit the result- 
ing model photometry to the data and compute the \ 2 statistic. Af- 
ter marginalising over all parameters, we obtain a stellar mass and 
an estimate of its error. Higher S/N spectra are required to constrain 
the mass ratio between with older and younger stellar components. 
Therefore our systematic errors are much larger than the computed 
formal statistical errors. To quantify this uncertainty, w e compared 



our m ass estimates using a second SED-fitting code l Kriek etlrtl 
120091) . again using a Salpeter IMF and lMarastor]dl998l.l2005h mod- 
els. The second code fits a single model template and therefore 
measures only a young stellar population component, excluding the 
additional free paramete r of a second older stellar population from 
ISchawinski et alj d2007l) . The single stellar population fits result in 
a total stellar mass 0.75 dex lower on average. This result is not 
unexpected as a younger stellar population is more luminous and 
can account for the same amount of light with a smaller mass con- 
tribution. Accounting for the average 0.75 offset between the two 
mass measurements, the residual dispersion was just under 0.3 dex. 
Therefore each individual mass measurement is uncertain at this 
minimum level. Because two stellar populations are more likely 
to reside in these galaxies (§ 15. U , we quote the masses from the 
two-burst model for the results in this paper acknowledging their 
uncertainty. Figure [10] shows an example stellar population fit. The 
median stellar mass of a Pea in our sample ~ 1O 9 ' 5 M0 and they 
range from 1O 8 ' 5 M0 to nearly 1O 1O ' 5 M0, indicating the Peas as a 
class are significantly less massive than an L* galaxy. Galaxy mass 
estimates are also included in Table[4] Column 9. 



4.3 Metallicity 

We measure gas phase metallicity for the Peas using [N II] 
A6584 A/ [O II] AA37 26, 3729 Ato estimate log[0/H] + 12 
dKewlev & Dopital l2002h . The average Peas has a metallicity of 
log[0/H] + 12 ~ 8.7. These metallicity measurements are broadly 
consistent with those determined from other line ratios for which 
we have lower (S/N) and fewer galaxies with measurements as well 
as with the location of the Peas on the BPT plot (Figure [6}. We 
plot th e Peas on the mass metallicity r elation fromlTremonti et all 
(2004) in Figure[TT] The lines are from lTremonti et all d2004l) . en- 
closing 68 % (dotted) and 95% (dashed) of the star-forming galax- 
ies in the lTremonti et al.l d2004l) sample. Although we do not see a 
trend in the Peas' metallicity with mass, they are roughly consistent 
with the mass-metallicity relation. The exception to this agreement 
is for the Peas with the largest masses, which have the same low 
metallicity as their lower-mass counterparts and therefore lie below 
the mass-metallicity relation. This is likely due to the uncertainty in 
the mass determinations rather than the measurement of gas phase 
metallicity. We include metallicity in Table[4] Column 8. 

Overall, we find the Peas have log[ Q/H] + 12 ~ 8.7, sub- 
solar (Zq ~0.5) as per measure ments of | Grevesse & Sauvall d 19981) 
and helioseismic measurements (Basu & Antiall2008l) . but near the 
solar abundances of Aspl und et alj d2005l) . These metallicities are 
common in low mass galaxies like the Peas. 



4.4 Specific Star Formation Rate 

Combining the mass measurements with the star formation rates, 
we find extraordinarily high specific star formation rates for the 
Peas. "Specific" star formation rate (SSFR) refers to the star forma- 
tion rate per solar mass in units of year - 1 and thus can be directly 
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Figure 11. The mass metallicity relation after lTremonti et al] 120041) . The 
Peas are sub-solar metallicity. Although the Peas are in general consistent 
with the mass-metallicity relation, they depart from it at the highest mass 
end and thus do not follow the trend. The Peas selection selects galaxies 
with a range of masses, but a more uniform metallicity. 



related to the time taken to double the stellar mass of a system 
(1/SSFR). In Figure [T2l we comp are the Peas to the Galaxy Zoo 
Merger Sample (Darg etaT]|200a) . Major mergers are frequently 
sites of active star formation, and yet the Peas are an order of mag- 
nitude higher in specific star formation rates than this comparison 
sample. These specific star formation rates imply doubling times 
between 100 Myr and ~ 1 Gyr. The uniformly high star forming 
rates of the Peas are not unexpected because their selection cri- 
teria targets strong emission lines. If the rates of star formation 
were lower, they would not be detected as "green" in the SDSS 
imaging. However, the SSFR we measure are unusually high for 
alaxies at z^0.2, which typically reach SS FR~ 10 _9 yr _1 at most 



Brinchmann et al . 2004; Bauer et al. 2005). 



4.5 UV Luminosity 

The Peas are well spread throughout the 8,400 square degrees cov- 
ered by the SDSS and not concentrated in any of the areas covered 
by deep multi-wavelength surveys. However, GALEX GR4 is well 
matched to SDSS in depth a nd area and 139 of the 251 Peas are 
detected in GALEX GR4 data dMorrissev et alj2007t) . For the 57 of 
80 star-forming Peas with Galex detections (S/N > 3), the median 
luminosity is 3 x 10 10 Lq. We include Lpuv in Table [4] Column 
6. The very high UV luminosities combined with low reddening 
(ED are rare in local galaxies and are more typically found in 
galaxies at higher redshift. 



5 COMPARISON WITH OTHER SAMPLES OF 
GALAXIES 

5.1 Blue Compact Dwarfs 

Blue compact dwarfs are a subset of the local dwarf pop- 
ulation, first identified by IZwickvl d 19651) as star-like field 
galaxies on Palomar Sky Survey Plates. They are character- 
ized by compact, gas-rich regions of high star formation rates 
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Figure 10. An example SED fit to a typical starforming Pea. The top panel shows both the original emission-line dominated SDSS ugriz photometry (large 
black points), and the medium-band photometry derived from the emission-line subtracted spectrum used for the SED fit (open boxes). The best-fit template 
(small black points and dashed line) is plotted on top. In the bottom panel, we show the SDSS spectrum of this Pea together with the ugriz filter transmission 
curves. The [OIII] line dominates the r-band flux and the i-band flux is significantly affected by Ha emission. The SDSS broad-band photometry is dominated 
by emission lines to such an extent that fitting stellar templates to it cannot yield reliable results. 
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Figure 12. Specific star-formation rate vs. mass for the Peas (purple dia- 
monds) and the Galaxy Zoo Mergers sample (black points). The Peas have 
low masses typical of dwarf galaxies and much higher specific star for- 
mation rates compared to the merging galaxies. The dashed line shows a 
constant star formation rate of 10 MQyr -1 . Most of the Peas had SFRs 
between 3 and 30 MQyr -1 and hence follow this line closely. 



(Papa deros et al.l l2002h IVaduvescu et alj 120071) . often lying in- 
side an older stellar population of larger spatial extent contain- 
ing a significant mas s contribution of the galaxy as a whole 
dLoose & Thuanl Il98rj ; Ikunth et al.l [l988l ; IPapaderos et~ai1ll996t 
lAloisi et al.l 120071) . They are generally low-metallicity (7.12 ^ 
12+log(0/H) ^ 8.4), evolved, gas-rich dwarfs undergoing recur- 
rent starburst activity i Loose & Thuanl 198fJ:JPapaderos et al .1 19961: 



iGil de Paz & MadoreH2005l ; IPapaderos et al.l2008l) . As a class Blue 
Compact Dwarfs cover a wide range in absolute magnitude and are 
often sub-classified by size (Ultra Compact Blue Dwarfs) or lumi- 
nosity (Luminous Blue Compact Galaxies). 



Ultra Compact Blue Dwarf galaxies are the smallest of the 
Blue Compact Dwarfs with physical diameters less than 1 kpc 
and low masses (~ 10 7 - 10 9 M^; ICorbin et al.l 120061) . Typi- 
cally an order of magnitude less luminous than other BCGs, their 
optical spectra can be dominated by very strong emission lines 
dGuseva et al.ll2004uCorbin et al J l 2006) and they can contain sub- 
stantial amounts of i nternal dust (eg. E(B — V) ~ 0.28 ; 
ICorbin & Vaccall2002l) . Like other low mass galaxies, Ultra Blue 
Compact Dwarfs have very low metal l icities , yet they still lie on 
an extrapol ation of the|Tremonti et ail d2004l) mass-metallicity re- 
lation dCorbin et alj|200q) . These galaxies can show asymmetric 
morph ologies with multiple sub-clumps of activ e star-forming re- 
gions dCorbin & Vaccall2002l;lc"orbin et alj|2006l) . While their light 
is dominated by the compact, young (~1-10 Myr) stellar popula- 
tion, their stellar mass is dominated by an older, evolved (~ 10 Gyr) 
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popul ation. As a popul ation, these galaxies tend to reside within 
voids dCorbin et all2006h . 

At the other end of the Blue Compact Dwarf category, Lumi- 
nous Blue Compact Galaxies are the most luminous members with 
Mb —17.5. They are more massive than their Ultra Compact 
count erparts, with average masses near 5x 10 9 Mq (Guzm an et al.l 
120031) . As one might expect from a class of more massive galax- 
ies, their metallicities are on average higher than the Ultra Blue 
Compact Dwarfs (7.7 s£ 12+log(0/H) sC 8.4; iGuseva et alj|2004 
IHqvos et alj|2007l) . However, like their Ultra Compact Blue Dwarf 
counterparts, Luminous Blue Compact Galaxies show compact 
prominent regions of star formation, often with disturbed mor- 
phologies iBerg vali & Johanssonl 1 19851) and star formation rates 
typically ranging from 1-5 M0yr _1 . Underlying these young, 
star-forming regions, older stellar populations with the colours 
and stellar profile s of older, massive ellip ticals are detected in 
near-IR imaging dBergvall & Ostlir] [2002). Examples of Lumi- 
nous Blue Compact Galaxies are found at intermediate re dshifts 
dGuzman et afl 1 19971 ; IPhillips et alj[l997l ; iHovos et alj|2004h and 
may contribute u p to 50% of the star formation rate density in the 
Universe at z=l dGuzman et aljl997h . 

Because both the Peas and Blue Compact Dwarf galaxies con- 
tain strong emission lines originating in compact regions of star for- 
mation, we investigated their potential overlap as a class. Although 
the Peas are similar in morphology, environment and physical size 
to the z=0 Ultra Blue Compact Dwarfs, they appear to be a different 
class of galaxies. The Peas have significantly higher metallicities, 
typically more than 0.5 Z© compared to 0.02 Zq for the Ultra Com- 
pact Blue Dwarfs. Additionally, the stellar masses of the Peas are 
on average ~10 9,5 Mq, roughly an order of magnitud e larger than 
simil ar measurements of Ultra Compact Blue Dwarfs dCorbin et al.l 
2006). While Ultra Compact Blue Dwarfs are found near z ~ 0, 
the Peas are detected at 0.1 13^z^0.36, meaning they have several 
Gyrs to grow and further increase their metallicity before reaching 
z=0, increasing the disparity in the measured masses and metallic- 
ities. In contrast, the largest of the Blue Compact Dwarfs, Lumi- 
nous Blue Compact Galaxies have stellar masses and metallicities 
that match the measured values for the Peas. The luminosities of 
the Peas are easily as bright as this class of galaxies, with median 
Mb ~ —20. Additionally, studies of low-metallicity SDSS objects 
like Blue Compact Dwarfs find that t hey lie in the sa me location 
on the BPT plot shown in Figure [6] (Izotov et al. 2006) and [O III] 
5007 / H/3 line ratios for a sample of four Luminous Compact Blue 
Galaxies measured with STIS (log([O777]A5007A 77/3) ~ 0.5; 
lHovosetal.ll2004l) are similar to those we find here (Figure [6J. 
Studies of the Peas at near-infrared wavelengths could potentially 
reveal older stellar populations, like those found in Luminous Blue 
Compact Galaxies. These older stellar populations could reveal ad- 
ditional stellar mass and potentially larger radii. Therefore, the Peas 
could potentially be classified as part of the heterogeneous Lumi- 
nous Blue Compact Galaxies category. 

We tentatively conclude that the Peas form a different class of 
galaxies than Ultra Blue Compact Dwarfs, but they are similar to 
the most luminous members of the Blue Compact Dwarfs category, 
Luminous Blue Compact Galaxies. 

5.2 Local UV-Luminous Galaxies 

iHeckman et al.l d2005l) defined a sample of 215 ultraviolet- 
luminous galaxies (UVLGs) from the SDSS spectroscopic galaxy 
sample (DR3) and GALEX (GR1) with luminosities Lfuv ^ 
2 x 1O 1O L . We have 139 Peas with Galex detections (S/N ^ 
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Figure 13. Neither the distribution of [O III] equivalent width for ~ 
10, 0000 comparison galaxies ( i)2.1l red lined histogram) nor the UVLG 
sample (blue solid histogram) have equivalent widths comparable to the 
Peas (green lined histogram). As a class the Peas have much larger [O 
III] equivalent widths than normal or UV-luminous galaxies. All equivalent 
widths displayed here are from the SDSS pipeline measurements. 



3), 44 of which meet Heckman et al. ( 2005)'s criteria (Lfuv ^ 
2 x 10 10 7/q, z<0.3 and SDSS pipeline spectral type identified 
as galaxy. Since GR4+DR7 cover a much larger area than 
GR1+DR3, we find onl y four overlapping s ources between the Peas 
and the UVLG sample dHoopes et al.l20070 . However, we note that 
none of the other 211 UVLGs fall into the SDSS colour-selection 
wedge of the Peas. 

The Peas selection includes many UV-luminous galaxies, but 
does not uncover the same population of galaxies as a selection 
based on UV luminosity. To illustrate this in Figure[T3] we compare 
the [O III] equivalent width, as measured by the SD SS pipeline, for 
the ga laxy sample described in Section l2~T1 and the Heckm an et al] 
(2005) UVLG sample to that of the Peas. Unsurprisingly, due to 
their selection, the Peas have much larger equivalent width mea- 
surements than either the UVLG sample or the general spectro- 
scopic galaxy sample (§[T3j. The Pea [O III] equivalent widths are 
due to the combination of both very high star formation rates and 
the faint continuum. The Peas form a sample of objects selected to 
have large [O III] equivalent widths, many of which also have large 
UV luminosities. As the UVLGs are currently one of the best local 
sa mples of galaxies ana logous to high redshift Lyman-break galax- 
ies dHoopes et al.l2007l) . we now look into the comparison between 
the Peas and high redshift galaxy samples. 

5.3 UV-luminous High Redshift Galaxies 

Peas are sites of extreme star-formation in the local Universe. They 
are both very compact and have low stellar masses, yet they have 
enormous star formation rates, as shown by their emission lines 
and their UV luminosity. This sample of galaxies is quite distinct 
from typical z ~ star-forming galaxies, i.e., blue spiral galax- 
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ies or dusty irregulars dKeimicutd[l998l) . Therefore, we look at the 
properties of higher redshift galaxie s where the bulk of t he star for- 
mation in the Universe is occ urring dMadau et al Jl 1 996l : iLill v et al.l 
ll99d : lReddv & Steidelll2009h . 

Various colour-selection techniques have been developed for 
selecting samples of high redshift galaxies, including the Lyman- 
break dropout technique, which selects z ~ 3 galaxies by identi - 
fying (7-band dropouts in deep imaging (e.g. JSteidel et al.ll2003h . 
and Lya emission selection, which uses narrow-band imaging to 
select galaxies at various redshifts due to the excess emission of 
Lyq compared to the neighbouring continuum ([C owie & Hu 1998; 
iThommes et"ai1 1 19981 : iKudritzki et all l200d : ISteidel et all l2000h . 
Both techniques return a similar population of galaxies, but Lya 
emitters are characterised by a fainter c ontinuum and st ronger 
equivalent widths in their emission lines jGiavaliscdl2002h . Al- 
though both selection techniques have different biases, the galaxies 
showi ng Lyq in emission t end to be smaller and younger on av- 
erage dGawiser et al. I l200l [Finkelstei n'etal]|2007h . UV-luminous 
high redshift galaxies, like Ly-a emitters and Lyman-break galax- 
ies, are characterized by a high UV luminosity and a relatively low 
obscuration by dust dGawiser et alj|2006l ; IVenemans et al.l 12005). 
They are starburst galax ies with strong blue continua dominated 
by young, massive stars jGiavalisccJ '2002). These techniques have 
been extended to other redshifts out to z ~ 6, u sing different pho- 
tometric bands to select drop-out galaxies (e.g.,|_O uchi et alj|200lh 
IStanwav et al.l2003l : lGlavalisco et al]2004l : lKakazu et al.l2007l) . In 
a broad sense, the Peas are an application of a Lya-like selection 
technique at much lower redshift and targeting the [O III] emission 
line. 

UV-luminous, high-redshift galaxies are similar in size to 
the Peas, with half light radii <; 2 kpc dGiavaliscoetalJll996l ; 
iBremer et alj2004|Pascarelle et alj 1 998). Ma ny of these high red- 
shift galaxies show di sturbed morphologies JRavindranath et aU 



120061 : Iloiz et alj 120061) , like those seen in the HST images of 
the star-forming Peas in Figure [7] In Lyman-break galaxies, rest- 
fra me equivalent widt hs of [O III] can reach up to hundreds of 
A dPettini et al.ll200ll) . although most are much lower. For Lya 
emitters, most are selected to have EW of Lyq ^ 20 A and can 
range up to ~240 A or more dGawiser et al. 200t3; lGronwall et"al] 
2007). In this way, the Peas are more like Lya emitters, selected 
for their strong emission lines, and the subset of Lyman-break 
galaxies with large emission lines. Reddening measures of Lyman- 
break galaxies are low, similar to those we measure fo r the Peas, 
E(B - V) < 0.2 dShaplev et al.ll2003l : lGiavaliscdl2002h . The Lya 
line is also easily s uppressed by dust, making Lya emitters sim- 
ilarly low in dust dGawiser etai] |2006). Additionally they have 
high star formation rates of a few to tens of solar masses per 
year dBarmbv et alj|2004l: iLehmer et al.ll2005l; ICoppin et al ll2007l ; 
ICarilli et al.l 120081 ; Ilvlannucci et alj 120091 ; IPentericci et alj 12009b . 
again similar to what we measure for the Peas. Peas are similar 
to UV-luminous high redshift galaxies in size, morphology, large 
emission lines, reddening and star formation rate. 

Masses measured for high redshift UV-luminous galaxies 
are slightly larger than those fo r the average Pea, ~ 10 9 Mp 







lQ n M iBarmbv et al.1 [20041; ICoppin et alj |2007|; lYabe et al.l 



2009; Mannucci et al. 2009; Penteri cci et alj l2009ri . In terms of 
SSFR vs mass, Lyman Alpha Emitters li e just below t he 10 
Mpyr" 1 line (Figure 12, dashed line) dGawiser et al] 120061 : 



ICastro Ceron et al.N2006h . While LBGs range from 10 M yr" 
to 100 M©yr _1 , at m asses below 10 9 .5 they are preferen- 
tially near lOOMayr" 1 dCastro Ceron etall l200d IShaplev et ail 
|2001|; iBarmbv et al. 2004). Metallicities of Lyman-break galax- 



ies are typically meas ured near 10-50% solar dPettini et al.ll200ll : 
iMannucci et al. 1 20091) . significantly lower than that of the lower 
redshift Peas. This is not surprising since the Peas have had sig- 
nificantly longer to enrich their gas. Finally, Lyman-break galaxies 
are kn own to be stron gly clustered and found in the densest re- 
gions dGiavaliscdl2002h in contrast to the Peas. The peas are found 
in lower density regions, have lower masses and smaller specific 
star formation rates than the UV-luminous galaxies found at high 
redshift. The smaller mass and lower density environment of the 
Peas is consistent w ith a picture of downsizing dCowie et al.ll 19961 : 
iThomas et al ]|2005h where smaller present day galaxies form their 
stars at later times in lower density environments. 

Understanding the evolution of starburst galaxies over cosmic 
time is central to understanding the build up of stars in galaxies. At 
high redshift (1.9 ^ z ^ 3.4), UV-luminous galaxies are responsi- 
ble for the f ormation of a large frac tion (~ 40%) of the present day 
stellar mass ( Reddy & Steidel 2009). Therefore, the Peas are poten- 
tially the remnants of a mode of star formation that was common 
in the early Universe. If that is the case, then the Peas are an ideal 
laboratory for understanding this mode of star formation, as their 
continuum properties are easily accessible with large ground-based 
telescopes. Additionally, their low redshifts allow optical and near- 
infrared telescopes to investigate potential underlying older stellar 
mass components. High redshift galaxies are both small and faint, 
even at HST resolution, but the galaxi es at z ~ 0.1 — 0.3 can be 
imaged at higher physical resolution dOverzier et al. 2008). Stud- 
ies of the morphologies of low dust, high star-forming galaxies 
in the local Universe can lend insight to the processes occurring 
at higher redshift where the morphologies cannot be as finely re- 
solved. The X-ray luminosites of known Lyman-break galaxies are 
~ 10 41 ergs/s, accessible a t z~0.2 with snapshots from Chandra 
dHorns chemeier et al. 2008), allowing for studies of their star for- 
mation rates at X-ray wavelengths. The Peas may be the last rem- 
nants of a mode of star formation common in the early Universe, 
and therefore an excellent laboratory for understanding that mode. 



6 SUMMARY 

We investigated a class of galaxies, known as Peas, discovered by 
the Galaxy Zoo project. These galaxies are characterised by a dis- 
tinctly green colour in gri imaging arising from a very large [O III] 
equivalent width. 

• 251 Peas were collected from the SDSS spectroscopic 
database based on a colour selection in the redshift range 0.112 ^ 
z < 0.360. 

• The Peas are unresolved in SDSS imaging, placing an upper 
limit on their physical radius of approximately 5 kpc. 

• The median environmental density around the Peas is less than 
two-thirds of that around normal galaxies. 

• The BPT spectral line diagnostic reveals that the majority of 
the Peas are star-forming galaxies, some of which show patchy 
morphology in HST imaging. 

• We uncover 8 new Narrow Line Seyfert Is from the SDSS 
archive. They lie below the M-cr relation, similar to other samples 
of NLSls. 

• From a sample of 80 star-forming galaxies with high S/N 
spectral measurements, we find that the Peas have very large star 
formation rates (up to 30 Moyr -1 ), low stellar mass (~ 10 9 ' 5 
M ), low metallicity (log[0/H]+12~8.7) and large UV luminosi- 



ties (~ 3 x 1O 1U L ). 
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• The Peas form a different class of galaxies than Ultra Blue 
Compact Dwarfs, but may be similar to the most luminous mem- 
bers of the Blue Compact Dwarfs category. Luminous Blue Com- 
pact Galaxies are similar to the Peas in their masses, morphologies, 
metallicities, luminosities and redshifts. It would be interesting to 
study the Peas at NIR wavelengths to see if they have older underly- 
ing stellar populations like those found in Luminous Blue Compact 
Galaxies. 

• The Peas share properties similar to local UV-selected sam- 
ples in Sloan, but uncover a different population with more extreme 
equivalent widths of [O III] emission line. 

• The Peas are similar to UV-luminous high redshift galaxies 
such as Lyman-break Galaxies and Lya emitters. However, these 
high redshift galaxies are higher in mass, lower in metallicity and 
found in the densest regions. The smaller mass and lower density 
environment of the Peas is consistent with a picture of downsizing, 
where smaller present day galaxies form their stars at later times in 
lower density environments. If the underlying processes occurring 
in the Peas is similar to that found in the UV-luminous high red- 
shift galaxies, the Peas may be the last remnants of a mode of star 
formation common in the early Universe. 
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Table 4. Properties of Pea Star-Forming Galaxies 



SDSSID a 


RA a 


Dec a 


z a 


[0 III] EW a 


r b 


SFR C 


log(0/H)+12 d 


Stellar Mass c 




(J2000) 


(J2000) 




A 


[10 44 ergss" 1 ] 


[Mgyr- 1 ] 




[Mq] 



587725073921409255 
588848899919446344 
587725576962244831 
587731187273892048 
587731513693503653 
587724233716596882 
587727179006148758 
587724241767825591 
587724240158589061 
587726032778559604 
587726032253419628 
588010360138367359 
587726102030451047 
587729155743875234 
587728919520608387 
587729229297090692 
587725818034913419 
587730774416883967 
587730774965354630 
587728906099687546 
587725550133444775 
588009371762098262 
588011122502336742 
588011103712706632 
588013384341913605 
587732134315425958 
587729777439801619 
587729777446945029 
587732152555864324 
587732578845786234 
587733080270569500 
588297864714387604 
587735695911747673 
587735696987717870 
587733441055359356 
588017605211390138 
588017114517536797 
588017116132540589 
588018090541842668 
588018090013098618 
588016878295515268 
587735661007863875 
588016892783820948 
587735663159738526 
588018055114784812 
588018055652769997 
588017570848768137 
587736915687964980 
587736915687375248 
587738410863493299 
587735349111947338 
587738570859413642 
587736940372361382 
587739153352229578 
587738947196944678 
587738371672178952 
588017978880950451 
587739408388980778 
588017977277874181 
587739406242742472 
587739828742389914 
587739652107600089 
587739721387409964 
587741600420003946 
587741421099286852 
587741532770074773 
587741817851084830 
587741391573287017 
587741392649781464 
587739648351076573 
587741490367889543 
587741532781215844 



146.242618 


-0.762639 


0.3002 


312.8 


1.04 


(0.21) 


12.45 


(0.52) 


8.78 


8 


195.546460 


-0.087897 


0.2255 


834.4 


0.19 


(0.04) 


4.47 


(0.22) 


8.75 


| 


261.776373 


59.817273 


0.3472 


651.5 




(...) 


24.98 


(1.48) 


8.71 


| 


351.413453 


0.752012 


0.2770 


319.6 


0.99 


( 0.06) 


11.11 


(0.85) 


8.70 




50.687082 


0.745111 


0.3043 


282.1 


1 .09 


( 0.27) 


14.22 


(0.80) 


8.83 


| 


22.292299 


14.992956 


0.2800 


354.0 


1.53 


(0.21) 


13.37 


(0.72) 


8.77 


1 


45.839226 


-7.989791 


0.1650 


886.0 




(...) 


8.72 


(0.64) 


8.7] 


| 


51.556792 


-6.586816 


0.1621 


822.9 




(...) 


11.44 


(0.41) 


8.76 


1 


54.949128 


-7.428132 


0.2608 


397.6 


2.00 


(0.14) 


28.96 


(1.68) 


8.78 




164.319700 


2.535293 


0.3028 


348.5 


1.35 


(0.14) 


8.75 


(0.57) 


8.76 


1 


191.097382 


2.261231 


0.2395 


1151.4 




(...) 


25.24 


(1.30) 


8.70 


1 


130.570630 


3.635203 


0.2194 


567.5 




(...) 


6.60 


(0.21) 


8.55 


^ 


236.787938 


3.603914 


0.2314 


891.8 




(...) 


9.83 


(0.61) 


8.74 


1 


173.265848 


65.228162 


0.2414 


475.2 


1.61 


(0.18) 


5.43 


(0.24) 


8.66 


1 


212.938906 


62.653138 


0.2301 


529.0 


1.45 


(0.16) 


12.80 


(0.47) 


8.67 


1 


234.405309 


58.794575 


0.2143 


851.6 


0.44 


( 0.06) 


6.32 


(0.20) 


8.58 


s 


235.209139 


57.411652 


0.2944 


173.3 


3.50 


( 0.24) 


16.81 


(1.18) 


8.91 


s 


339.396081 


13.613062 


0.2936 


438.3 




(...) 


24.78 


(0.86) 


8.55 


| 


6.716985 


15.460460 


0.2136 


754.2 




(...) 


3.32 


(0.21) 


8.81 




117.403215 


33.621219 


0.2733 


339.2 


2.04 


( 0.26) 


58.83 


(2.61) 


8.79 


| 


156.563375 


63.552363 


0.3338 


933.3 


1.27 


( 0.28) 


3.76 


(0.28) 






170.582224 


61.912629 


0.2045 


960.7 


0.40 


( 0.07) 


5.59 


(0.23) 


8.70 


8 


181.772142 


61.586621 


0.2620 


760.2 


0.68 


(0.13) 


13.02 


(0.62) 






226.616617 


56.450741 


0.2786 


341.9 


1.08 


(0.17) 


23.88 


(0.93) 


8.90 


8 


141.501678 


44.460044 


0.1807 


651.5 


1.55 


(0.11) 


14.35 


(0.89) 


8.62 


i 


195.368010 


51.080893 


0.3479 


554.4 


3.20 


( 0.26) 


23.17 


(1.74) 


8.63 


1 


204.299529 


-2.434842 


0.2737 


399.2 




(...) 


8.63 


(0.35) 


8.76 


1 


220.630713 


-2.164466 


0.2938 


1456.6 


1.61 


(0.19) 


20.22 


(0.80) 


8.57 


1 


116.991682 


23.609113 


0.1552 


874.9 




(...) 


3.26 


(0.07) 






157.912214 


7.265701 


0.2525 


763.9 


0.57 


(0.18) 


6.04 


(0.20) 


8.52 


8 


163.378431 


52.631353 


0.2526 


418.2 


3.61 


( 0.26) 


27.59 


(1.09) 


8.78 


1 


131.975356 


33.615227 


0.3063 


323.7 


1.89 


( 0.20) 


21.46 


(0.71) 


8.81 


1 


204.919632 


55.461137 


0.2291 


42.1 


1.41 


(0.16) 


3.22 


(0.13) 


8.7] 


1 


213.630037 


54.515587 


0.2270 


773.7 




(...) 


4.04 


(0.10) 


8.58 


1 


251.527242 


31.514859 


0.2907 


868.4 


1.20 


(0.14) 


6.56 


(0.36) 






154.513517 


41.105860 


0.2371 


1191.4 


0.88 


(0.14) 


8.83 


(0.25) 


8.68 


8 


216.023868 


42.279524 


0.1848 


1348.8 


1.49 


(0.12) 


19.66 


(1.09) 


8.78 


| 


228.535985 


38.868716 


0.3324 


632.5 


1.77 


( 0.22) 


6.02 


(0.36) 






235.755108 


34.767079 


0.1875 


673.9 




(...) 


2.02 


(0.05) 






251.898063 


22.783002 


0.3138 


578.8 




(...) 


4.20 


(0.49) 


8.61 


8 


137.879799 


31.457439 


0.2718 


426.3 


1.05 


(0.12) 


9.18 


(0.37) 






139.260529 


31.872384 


0.3002 


219.7 


2.08 


(0.18) 


20.29 


(1.51) 


8.95 


8 

K 


148.712329 


37.365500 


0.2834 


279.9 


0.91) 


(0.16) 


8.42 


(0.50) 


8.78 


§ 


149.415718 


37.702114 


0.2867 


235.3 


1.96 


( 0.23) 


15.02 


(0.63) 






220.041419 


46.326930 


0.3008 


304.4 


3.17 


( 0.26) 


31.50 


(2.70) 


8.72 


8 


223.648271 


45.482288 


0.2687 


463.0 


1.23 


(0.13) 


21.17 


(2.67) 


8.84 


S 


192.144310 


12.567480 


0.2634 


1055.8 


1.30 


( 0.04) 


15.01 


(0.75) 


8.53 


H 


241.152768 


8.333082 


0.3123 


1388.3 


0.55 


(0.13) 


10.64 


(0.58) 


8.75 


§ 


239.858241 


8.688655 


0.2970 


904.8 


0.66 


(0.13) 


3.44 


(0.38) 


8.64 


1 


152.987850 


13.139471 


0.1439 


2388.3 




(...) 


7.99 


(0.22) 


8.64 


§ 


184.766599 


15.435698 


0.1957 


1487.9 




(...) 


7.47 


(4.83) 


8.82 


| 


204.867933 


15.278369 


0.1921 


1289.0 


0.64 


( 0.06) 


18.81 


(1.57) 


8.83 


8 


217.614622 


34.154720 


0.1911 


586.8 


0.31 


( 0.06) 


3.92 


(0.15) 






117.990764 


16.637010 


0.2647 


440.1 


0.66 


(0.13) 


6.29 


(0.21) 


8.65 


8 

H 


123.966679 


21.939902 


0.1410 


1582.0 


0.27 


( 0.03) 


3.96 


(0.05) 


8.55 


I 

8 


125.698590 


22.695578 


0.2163 


1040.6 


0.90 


(0.13) 


37.41 


(4.15) 


8.81 


8 

H 


150.556494 


34.704908 


0.3210 


222.1 




(...) 


4.43 


(0.26) 


8.64 


§ 


174.342249 


35.407413 


0.1945 


647.2 


1.88 


( 0.08) 


20.38 


(0.91) 


8.76 


8 


171.657352 


38.050810 


0.2469 


646.9 


1.56 


(0.11) 


24.18 


(0.64) 


8.56 


§ 
8 


178.020352 


34.013853 


0.3420 


1095.0 




(...) 


15.61 


(0.91) 






224.396405 


22.533833 


0.1488 


1563.1 


0.45 


( 0.06) 


8.94 


(0.29) 


8.55 


8 


238.041673 


21.053410 


0.2332 


893.5 




(...) 


7.29 


(0.44) 


8.68 


1 


249.330431 


14.651378 


0.2923 


452.5 




(...) 


9.89 


(0.17) 


8.62 




181.252807 


26.346595 


0.3427 


305.2 


1.06 


( 0.24) 


18.45 


(0.94) 


8.81 


§ 


126.715863 


18.347732 


0.2972 


774.4 


0.57 


(0.18) 


4.53 


(0.20) 


8.75 




133.350354 


19.506280 


0.2365 


860.3 




(...) 


7.50 


(0.41) 


8.68 




137.805603 


18.518936 


0.2622 


329.8 


1.72 


(0.16) 


26.17 


(1.67) 


8.89 




145.946756 


26.345161 


0.2366 


12.6 




(...) 


3.15 


(0.07) 


8.72 




152.329151 


29.272638 


0.2219 


1042.8 


0.34 


( 0.09) 


4.51 


(0.19) 


8.69 


1 


155.239428 


29.624017 


0.2555 


378.8 




(...) 


6.17 


(0.25) 


8.82 


5 


158.112322 


27.298680 


0.1924 


840.9 


0.56 


( 0.09) 


12.74 


(0.41) 


8.69 


1 

8 


159.779860 


27.472509 


0.2801 


182.2 


0.81 


(0.11) 


4.57 


(0.22) 
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Table 4 (cont'd) 



SDSSID a RA a Dcc a z a [0 III] EW a L Fuv b SFR c log(0/H)+12 d Stellar Mass" 

(J2000) (J2000) A [10 44 crgss" 1 ] [M yr _1 ] [Mq] 



587742014876745993 


141.869487 


17.671838 


0.2883 


760.7 


0.86 


(0.16) 


12.38 


(0.36) 


8.51 | 


68 
7<A 


9.38 


588023240745943289 


140.705287 


19.227629 


0.3175 


526.8 


0.45 


(0.13) 


14.38 


(0.60) 


8.59 1 


10.48 


587745243087372534 


141.384863 


14.053623 


0.3013 


1430.4 


0.97 


(0.13) 


18.73 


(1.13) 


8.73 | 


11 

% 
59 


9.08 


587742628534026489 


243.276317 


9.496990 


0.2993 


215.7 




(...) 


14.14 


(0.72) 


8.70 1 


9.65 


587744874785145599 


121.325174 


9.425978 


0.3304 


534.6 




(...) 


21.67 


(1.68) 






10.24 


587742013825941802 


197.653081 


21.804731 


0.2832 


273.8 


1.08 


(0.10) 


9.70 


(0.68) 


8.83 % 


90 
75 


9.27 


587742062151467120 


196.734804 


22.694003 


0.2741 


879.9 


0.38 


(0.12) 


8.83 


(0.32) 




10.10 


587741727655919734 


193.761316 


25.935911 


0.3119 


660.3 


0.99 


(0.13) 


7.20 


(0.43) 






10.02 



' l The values presented here are from the SDSS DR7 archive. 

b The values presented here are computed from the GALEX GR4 archive Far UV fluxes using the redshifts in Column 4. 

c The star formation rates presented here are from the Ha line. Errors are computed from 1 a flux errors. See § |4. 1| 

The metallicities shown here are calculated from the [NII]/[OII] ratio JCewley & Dopita 20021. The metallicity is followed by upper (above) and lower (below) 1 a errors. 
See §[43] 

c The masses are calculated following the methodology of Schawinski et al. ' 2007). See § 14.21 
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